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We study the instabilities of a harmonic oscillator subject to additive and dichotomous multiplica- 
tive noise, focussing on the dependance of the instability threshold on the mass. For multiplicative 
noise in the damping, the instability threshold is crossed as the mass is decreased, as long as the 
smaller damping is in fact negative. For multiplicative noise in the stiffness, the situation is more 
complicated and in fact the transition is reentrant for intermediate noise strength and damping. 
For multiplicative noise in the mass, the results depend on the implementation of the noise. One 
can take the velocity or the momentum to be conserved as the mass is changed. In these cases 
increasing the mass destabilizes the system. Alternatively, if the change in mass is caused by the 
accretion/loss of particles to the Brownian particle, these processes are asymmetric with momentum 
conserved upon accretion and velocity upon loss. In this case, there is no instability, as opposed to 
the other two implementations. We also study the distribution of the energy, finding a power-law 
cutoff at a value which increases with time. 
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I. INTRODUCTION 



The harmonic oscillator is, of course, one of the most basic models in physics, and the noisy harmonic oscillator is 
a touchstone of statistical physics. The simplest model of a noisy oscillator is that of an oscillator driven by additive, 
white noise: 



rn - 



2^ + kx = r](t) (1) 



dt 2 1 dt 

with white noise r\ (£), 

(»?(*)) =0, ( V (t 1 )r,(t 2 )) = S6(t 2 -t 1 ) (2) 

which, as is well known, leads to a stationary Gaussian distribution for x and v = dx/dt, whose width is proportional 
to \/~S. It is also possible to add noise, here multiplicative noise, via external modification of one of the system 
parameters. One such a model is that of a harmonic oscillator with random damping strength, j(t): 



m^f + 27(1 + + /,.,• =0 !•'-!) 



with another white noise. This equation was first used [T] to analyze water waves influenced by a turbulent 
wind field. It also can be used to describe, upon replacing x and t by the order parameter and spatial coordinate, 
respectively, phase transition dynamics in a moving system [2 . Here the first derivative term represents advection by 
some external flow. In this guise, the equation has been used to study phase transitions under shear [3], open flows 
of liquids [3] , dendritic growth [5] , chemical waves [Hj and motion of vortices [7] . 
Another possibility is to randomly modulate the stiffness parameter, k, so that 

m^+2^ + k(l + m)x = 0- (4) 

There are also many applications of this model in different fields in physics, such as wave propagation in a random 
medium [5] , spin precession in a random external field [9] , turbulent flow on the ocean surface |10j , and as well as in 
biology (population dynamics [H]), in economics (stock market prices [12]) and so on. 

Recently, one of us 13 - 16J considered yet another way of introducing multiplicative noise, namely via a fluctuating 
mass term. Here, as opposed to the other two possibilities outlined above, the mass term must always be positive, 
so white noise is ruled out. One possibility is dichotomous noise, where the mass alternates randomly between two 
values, m±. so that the equation reads 

m±^ + 2 7 ^ + fc(l + = 0. (5) 

This model was originally conceived off in the context of a Brownian particle undergoing random adsorption and 
desorption. Many other applications of an oscillator with a random mass have been considered |17j , including ion- ion 
reactions [H]-[2S], electrodeposition |5T], granular flow [S^'EI], cosmology [25]-[27]) film deposition [35], traffic jams 
[29]- [30], and the stock market [S]- [32]. 

The most striking phenomenon of these models with multiplicative noise is the existence of transitions, with no 
steady state for sufficiently large noise. It has been traditional in these studies of harmonic oscillators subject to 
multiplicative noise to scale out the mass, setting it to unity. This is of course possible, and one can write the critical 
noise amplitude (which is dimensionless) in terms of the dimensionless parameters k/(mX 2 ) and j/(m\), where A 
is the inverse correlation time of the noise. (We will discuss the parallel situation in the white noise limit below). 
However, doing so obscures an important aspect of the physics of the problem; namely, the nontrivial dependence on 
the mass. Thus, the mass is an important parameter in a way that it is not in the presence of additive noise only. 
One major goal of this work is to consider the mass dependence explicitly. We will see that in some cases the mass 
dependence is even nonmonotonic. In addition, we explore some other features of the problem which seem to have 
escaped notice. Primary among them is the fact that the statistics of the observables, x, v and E are anomalous, 
with power-law tails. As opposed to models with restabilization at a nonlinear fixed point [33 , here there is no cutoff 
of the power-law at long times. This is related also to the scale invariance of the harmonic oscillator. Thus, even for 
parameters for which (E), say, is finite, higher moments of E are infinite in the infinite-time limit. They are finite 
at finite time due to the presence of a time-dependent cutoff on the power-law distribution, so that the moments 
grew exponentially in time. This behavior is also reflected in a single time trace, with the dynamics having a bursty, 
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intermittent character |33) . The exponent governing the tail of the distribution determines which moments have 
finite infinite-time limits. Thus the power-law nature of the distribution is intimately related to the infinite series of 
transitions in the model, [34] corresponding to the onset of convergence of higher and higher moments. 

The plan of the paper is as follows. We first discuss the mass dependence of the well-known transition corresponding 
to the divergence of the average energy (equivalently (a; 2 ))- We first discuss the case of random damping, where most 
of the general features are already present. We then move on to treat the cases of random stiffness and random mass 
in turn. In the second section, we discuss the distribution of the energy and their connection to the infinite set of 
transitions for different moments. In the last section we present our conclusions. 



II. THE MASS DEPENDENCE OF THE ENERGY INSTABILITY 



A. Random Damping 



We first discuss the case of random damping. We work in the framework of symmetric dichotomous noise, recovering 
the white-noise case by taking the small correlation time, large amplitude limit. Thus the damping 7 switches randomly 
in a Poisson fashion between two values, j±. If both j± are positive, there is of course no instability. However, a 
negative 7_ by itself is not enough to ensure instability. The basic starting point of the analysis is the equations of 
motion of the second moments. These moment equation can be derived either from the Langevin equation or from 
a Fokker-Planck formalism. The former method was presented, for example, in Ref. [T3TU6| . In the Fokker-Planck 
formalism, we construct two coupled Fokker-Planck equations for the probability density P±(x,v), conditioned on the 
value of 7. The equations read 

D d 2 2 7+ d , . kx d d A . 

p + = —jr^ p + + — T~ vP + + — fT p + - v jt p + + 9 p + - p -) 

m z ov z m ov m ov ox I 

D d 2 2 7 _ d , N kx d d A 

ui ov z m Ov m ov ox 2 

Here, the diffusion constant D is related to the strength of the additive noise and A is the decay rate of the correlator 
of the dichotomous noise, and is inversely proportional to the average switching time, t, by A = 2/r. From these two 
equations, it is easy to generate a closed set of moment equations for the six conditioned quantities, (x 2 )±, {v 2 )± and 
(xv)±. The equations read 

\ ((x 2 ) + - <* 2 >_) 

2 A 

l±{xv)± + (v 2 ) ± T- - ((xv) + - (xv)-) 

m I 

xv) ± -- 1± (v 2 ) ±T U(v 2 ) + ~(v 2 )-) (7) 

These are, of course, equivalent to the moment equations arising from the Langevin equation. Setting the time 
derivatives to zero and solving yields 

= 1 2 v , ,_3v \ = D ( 7+ + 7-)(4fc + 3mA 2 ) + AXj^ + 7 2 ) + m 2 X 3 

y " ' 2 {y /+ y fc ( 7+ +7_)(m 2 A 3 +4mAfc + 8A7 + 7_) + (7 2 +7 2 )mA 2 +7+7_(10mA 2 + 16A:) K > 

Setting the mass to unity yields the result presented in earlier works. Alternatively, our present result can be 
reconstructed from the standard, m = 1 result by the substitutions 

7±^; D ^R_ (9) 

m m ui 

It is clear that as long as j± are both positive, (a; 2 ) is positive, and there is no instability. It is convenient to rewrite 
(x 2 ) in terms of 7 and a, where y± = 7(1 ± a), yielding 

2 D 87/c + 4fcmA + 6m7A 2 +m 2 A 3 + 87 2 A(l + cr 2 ) 

yX ' ~ 2k^ 8 7 fc + 4fcmA + 6m 7 A 2 + m 2 A 3 + 8 7 2 A - 47cr 2 (mA 2 + 2k + 2 7 A) ' ' 

Let us first consider the white noise limit: 

2a 2 

A — > 00 ; a 2 — > 00 ; — — = S fixed (11) 
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FIG. 1. Instability boundaries for the harmonic oscillator with random damping. Above the uppermost line, all moments of the 
energy diverge. Between this line and the curve underneath, the first moment of the energy is finite, and all higher moments 
diverge. Similarly, below the bottom line, the second moment of the energy is also finite. Below the bottom curve, the third 
moment of the energy is finite as well. Note that for a < 1, i.e. both y± > 0, there is no instability. Here, 7 = 0.025, k = 1, 
and A = 0.2. 



where (x ) reduces to 



D 



WN : 



fc(l - 2Sj/m) 



(12) 



Here the critical noise amplitude is S c = 771/(27), beyond which no steady-state exists. This critical noise amplitude 
is proportional to the mass, so that low mass systems are more unstable. 

Turning now to the general dichotomous system, we see that the numerator is positive definite, while the denomi- 
nator is negative for sufficiently large a 2 . The critical a 2 is 



2 87k + Akm\ + 6 7717 A 2 + to 2 A 3 + 
o r = — , „ ——, — 



4 7 (mA 2 + 2k + 2 7 A) 



8 7 2 A 



(13) 



It can be checked that cr 2 (m = 0) = 1 and increases monotonically with m, diverging for large m as m\/(4r/), in 
accord with our white noise limit result. This implies that for any given a 2 > 1, the system is unstable for small m 
and is stable for m > m c (cr), where the critical m, m c , increases with a. The instability boundary is indicated by 
the upper curve in Fig. [l] (the other curves will be discussed below in Sec. II). We thus see that for the dichotomous 
noise limit, the mass dependence is quite complicated. 

As already pointed out by Kitahara, et al. [35], an equivalent way to demonstrate the instability is to examine 
the dynamics of the system without external additive driving 77. Then, the system Eqs. Q is a homogenous linear 
system characterized by a 6 x 6 matrix. As usual, the presence of positive eigenvalues corresponds to existence of an 
instability, which yields precisely the same stability criterion as given above. 
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B. Random Spring Stiffness 



We now turn to the case where the spring stiffness, k, varies between two values, k±. The Fokker-Planck equation 
and the resulting moment equations are directly parallel to the previous case, and the steady-state solution is 

(r 2 )=4n (mA + 27) [2(fc + + fc_)+4 7 A + rnA 2 ] 

V ; (32 7 mA + 64 7 2 )A: + A;_ + (247 2 mA 2 + 47m 2 A 3 + 327 3 A)(fc + + fc_)-(fc + -fc_) 2 TO 2 A 2 1 ' 

As before, this result can be recovered from the standard m = 1 version by the substitutions given in Eq. ^. The 
first thing to note is that, as opposed to the previous case, here the denominator can be negative even if k± are both 
positive. This is due to the fact that changing for fe_ to k+ increases the energy of the system, and even though the 
reverse transition decreases the energy, the two do not have to balance. Rewriting things in terms of k and tr, with 
k± = k(l ± a), yields 

2 D (mA + 27) (4k + 4 7 A + mA 2 ) 

' ~ ~k 2 7 (mA + 2 7 ) (mA 2 + 4A 7 + 4k) - a 2 k(m\ + 4 7 ) 2 ' 15 -* 

In the white noise limit, this reduces to 

{x2)wN= 2 1 k(l- D kS/(4 1 )) (16) 

Thus in the white noise limit, the system is unstable for S > 4 7 /fc independent of the mass, as opposed to the case of 
the white noise random damping. 

Reverting to the general dichotomous case, we find that the critical <j 2 {m) above which the instability sets in is 

2= 2 7 (mA + 2 7 )(mA 2 +4A 7 + 4fc) 
° c fc(mA + 4 7 ) 2 1 ' 

In general, this function is not monotonic in m. It has an extremum at 

4 7 2 

m PT f = — (18) 

so that for 7 < 2k/ X, it has a maximum at positive m, with the value 



2 

a ext 



1+ ^r 

2k 



(19) 



In addition, note that 



a c 2 (0) = l + ^; 0^(00) = ^ (20) 



Thus, for 7 > 2k/X, a 2 is monotonically increasing in m, whereas for k/X < 7 < 2k/X, a 2 rises with m and then fails 
to a value greater than its value at 0, and for 7 < k/X, a 2 rises with m and then falls below its m = value. The 
upshot of this is that for large damping, 7 > 2k/ X, the system is stable for a 2 < er 2 (0), independent of mass, while 
for larger er 2 , the system is stable only above some critical value of m which increases with a 2 . For smaller damping, 
the situation is more complex. For intermediate damping, k/X < 7 < 2k/ X, the system is always stable independent 
of m for a 2 < a 2 (0). For cr 2 (0) < a 2 < a 2 (00), the system is stable above a critical mass which increases with a 2 . For 
a 2 (00) < a 2 < a 2 ^ the system is stable in a finite range of m, which shrinks from top and bottom as a 2 increases. 
Finally, for a 2 > o™*, the system is always unstable, independent of m. For 7 < k/X, on the other hand, the system 
is stable for a 2 < a 2 (oo), while in the very narrow band a 2 (00) < a 2 < <J 2 xt , the system is stable for an intermediate 
range of m's. For a 2 > cr 2 xt , the system is always unstable. The cases of very small and large damping are shown in 
Figs. [2J It should also be noted that the system is unstable with a 2 < 1, i.e. with both k± > 0, only for very weak 
damping, 7 < fc/(2A), and for sufficiently large m, 

m > 2 f 1) 37A+V7A(4 fc + 7A) 

\XJ k-2-yX K ' 
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FIG. 2. Instability boundaries for the harmonic oscillator with random stiffness. Above the uppermost line, all moments of the 
energy diverge. Between this line and the curve underneath, the first moment of the energy is finite, and all higher moments 
diverge. Similarly, below the bottom line, the second moment of the energy is also finite. Below the bottom curve, the third 
moment of the energy is finite as well. Note that for in the left panel, with 7 small, there is an instability at large m even for 
a < 1 , i.e. both k± > 0. In the left panel, 7 = 0.4. In the right panel, 7 = 3. In both cases, k = 1, and A = 1. 



C. Random Mass 



Recently )13H16j . a model in which the mass term randomly fluctuates was studied. In this case, it is clearly 
important that the fluctuations do not change the sign of the mass term, which rules out the possibility of white 
noise. Nevertheless, it is possible to consider the case of dichotomous noise where the "mass" fluctuates between two 
values, m±. Depending on the system being modeled, the random "mass" can represent, e.g., the inductance in an 
LRC circuit, the moment of inertia in a physical pendulum, etc. Depending on the physical application, different 
mathematical models are called for. For example, in the case of a random inductance, the charge and current are 
conserved at the moment of the fluctuation. For the case of a random moment of inertia, one would expected the 
angle and angular momentum to be conserved at the moment of the fluctuation. Thus, these two cases lead to two 
different mathematical models, one which reads 



d 2 

m± 

and the second which reads 



m± -tttX = — 7± — kx + i] (22) 
dt z 



— (m±x) = — 7± — kx + r\ (23) 

In the first, continuous velocity, case, which was the case studied in Refs. |13H16j . the moment equations read 

2^{x 2 ) ± =2{xv) ±T \{{x 2 ) + - {x 2 )S) 

~{xv)± = — —{x 2 ) ± - — (xv)± + (v 2 } ± T ~ ((xv} + - 
at m+ m+ 2 



The solution is 
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Uv 2 )± = ™ - ~{xv) ± - ^V) ± T £ «- 2 } + - <« 2 }-) (24) 
at m± m± m± 2 



2 2D \(k + 7A)(m + + m) 3 + 167(fc + 7A)m+m_ + r )\ 2 2m + m-(m + + m_) + A 3 m + m_(m^ + m 2 _) 

c ) = 
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(25) 
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FIG. 3. Instability boundaries for the harmonic oscillator with random mass coefficient, and velocity (or momentum, the phase 
boundaries coincide in the two cases). Above the uppermost line, all moments of the energy diverge. Between this line and the 
curve underneath, the first moment of the energy is finite, and all higher moments diverge. Similarly, below the bottom line, 
the second moment of the energy is also finite. Below the bottom curve, the third moment of the energy is finite as well. Here 
the model is physical only for a < 1 , i.e. both m± > 0, there is no instability. In the left panel, 7 = 0.1. In the right panel, 
7 = 0.3 Here, k = 1, and A = 1. 



As opposed to the cases of random damping and random stiffness considered above, here for physical reasons the 
mass must always be positive, and so we cannot pass to the white noise limit. Examining the dichotomous result, 



Eq. (25), the numerator is positive definite for m± > 0, while the denominator vanishes at a critical noise strength, 



a, where m± = mil ± a): 

2 2 7 (2 7 + mA)(4 7 A + 4fc + mA 2 ) 
" c = ^kX 2 (26) 

It is clear that a 2 is a monotonically decreasing function of m, approaching 2 7 A/fc as m gets large. Since the positivity 
of m± requires that a 2 < 1, there is clearly no instability if 2 7 A > k. If 2 7 A < k, then there is an instability above some 
critical value of m, if a 2 is large enough such that u 2 (m) < a 2 < 1. Turning this around, for fixed 1 > a 2 > 2 7 A/fc, 
there is an instability for m greater than a critical value. For m < 7 /A(vl + 16r + 32r 2 — 3 — 4r)/(l — r), where 
r = fc/(2 7 A) < 1, there is no instability for any a 2 < 1. This is illustrated in Fig. [3J 

For the second case, with momentum conserved at the transition, the moment equations read 





771 zp 


\ (m± 






m± 






(<»*>*-( 






, m ±) 



j t {x 2 ) ± =2{ X v) ±T \({x 2 ) + ~{ X 2 )^ 

|m ± = -A ( , 2)± _|l m± + (u2)± 

|(« 2 > ± = % - —{*v) ± - ^V) ± - } - [ (v 2 ) ± - ( — ) (v% I (27) 
at m± m± m± I \ \ m ± J J 

One way to see that these are the correct equations is to write the moment equations for the variables {(x 2 }± , (xp}±, (p )±} 
where p± = m±v±. In this case, the solution for (x 2 ) is 

2 D (4 7 + A(m + + 777_ ))(8fc + 8 7 A + A 2 ( 777+ + 777_)) 

[ ' ~ ~k U-/ 2 {k + 7 A) + (m+ + m_)(16/c 7 A + 24 7 2 A 2 ) + [m\ + m 2 )(2 7 A 3 - k\ 2 ) + rn + m_(4 7 A 3 + 2A 2 ) 1 ' 

As expected, the numerator is positive definite for m± > 0, while the denominator is negative if 

a 2 > 2 7 (2 7 + mA)(4A 7 + 4fc + mA 2 ) 
m 2 fcA 2 
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This is exactly the same condition for instability as for the conserved velocity model, even though (x 2 ) is different in 
the two models. The same is true for the instability conditions for the higher-order moments. 

There is yet a third possible model. If the changes in the "mass" term represent actual changes in the mass due 
to accretion/desorption of particles to the Brownian particle, then due to Newton's third law, the addition of mass 
conserves momentum, whereas the loss of mass conserves velocity. Then the moment equations read 

d k 2 27 2 A / , rn_ 

— (xv)+ = {x )+ (xv) + + (v )+ - - (xv)+ (xv) 

at m+ m+ 2 \ 

Jt {XV) - 

d 2 'ID -Ik , 4- , A , / , , 

di {v)+ = ^-^ {xv)+ -^ {v)+ -2 [ {v)+ -{^; ] 1 ' 
i^- = ^ - ^<->- - ttw- - 2 C<0- - <o + ) (30) 

In this model, the Brownian particle losses kinetic energy both upon accretion and loss and thus the pumping 
mechanism present in the other versions is absent. Indeed, solving the FPEs shows that the system does not possess 
an instability for any m± > 0! 
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III. THE DISTRIBUTION OF ENERGY AND HIGHER TRANSITIONS 



A. Random Damping 

Examining the how the distribution of E = (kx 2 + mv 2 )/2 changes with time is very instructive. We see in Fig. 
[4] that the distribution is characterized by a power-law tail which is cut off beyond some value of E, a value which 
increases with time. Given that in the absence of such a cutoff the mean value of E would diverge (the power is 
smaller than 2 in magnitude), the mean is strongly time-dependent, as the moment equations predict. A similar 
phenomenon happens for a random walk in a logarithmic potential |36j . where the distribution of position also has 
a time-dependent cutoff (that in that case grows as the square-root of time). As we lower the noise amplitude, or 
increase the mass, the power increases in magnitude. The critical noise amplitude, or mass, is that value for which 
the power is equal to -2, at which point the mean value of E is finite in the t — > 00 limit, when the power-law is 
not cut off. We see this in Fig. [5j where the power is close to -2 for m — 2, as this m is very close to the critical 
rn c = 2.045 for these parameters. This implies that just beyond the critical mass, the mean value of E converges but 
the mean value of E 2 diverges. This can be verified directly by looking at the moment equations for the fourth order 
moments, {(x 4 ) ± , (x 3 v)±, (x 2 v 2 )±, (xv 3 )±, (v 4 )±}. These moment equations read: 

|(x 4 ) ± =4( OT 3 ) ± T^((^ 4 ) + -(x 4 )_) 

^ (x 3„ )± = _ A (a; 4 )± _ ! 7±<a .3„ )± + 3{x 2 v 2 )± T A ( 3 _ 3 ^ 

at m m 2 

j t {x 2 v 2 ) ± = g - 2 A (X 3 U)± _ ± 7± (,V) ± + 2<x«'> ± T - 2 «*V) + - (.V)_) 

Uxv*) ± = ™{xv) ± - 3-(x 2 v 2 ) ± - 7± (^ 3 )± + (v 4 )± T $ ((x 3 v) + - (x\)„) 
at m z m m 2 

* {v * )± = ^(x 2 ) ± - ^(xv*) ± - ^ 7± (« 4 ) ± T £ ((. 4 ) + - (v 4 )-) (31) 
at m z m m 2 

These fourth-order moment equations are driven by the second-order moments. To locate the transition, it is preferable 
to look at the eigenvalues of the 10 x 10 stability operator. The critical value of m for which the fourth-order moments 
exists is to 4 = 4.132. Thus, we expect that the asymptotic slope of P(E) for m — 4.13 should be —3, as is also seen in 
Fig. [5] Thus, there is no steady-state solution of the fourth-order moment equations even for values of m for which 
a steady-state solution of the second-order moment system exists. Similarly, by looking at the sixth order moments, 
the condition for the existence of a steady-state solution is even more restrictive. This is shown in Fig. [TJ where 
the critical lines for the existence of stable solutions of the second, fourth and sixth moment equations are shown, 
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FIG. 4. The probability density of E for the harmonic oscillator with random damping for various times, as measured by 
simulation. Notice that as t increases, the large E behavior approaches a pure power-law, whereas for shorter times, it is 
cut-off. Here k = 1, A = 0.2, m = 1, 7 = 0.025 



corresponding to the borders of the regions of finiteness of (E), (E 2 ) and (E 3 ). The critical a for the nth moment 
starts off linearly from a = 1 for m = 0, turning over and behaving as 



2Ato f 32 ) 
3n + 2 {62) 

for large to. 

These anomalous statistics for E imply that the system has an intermittent, bursty character. We see this in Fig. 
6j where we plot E(t) for a single long run. The figure demonstrates clearly that the system is quiescent for a while, 
then undergoes after some random wait a large fluctuation, and returns to its quiescent state. 
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B. Random Spring Stiffness 

As with the case of random damping, here also the fourth order moments are unstable for a wider range of 
parameters than the second order moment system. Likewise, the sixth order moments are unstable for yet a wider 
range of parameters. This is shown in Fig. [2] The to = value of the critical a for the general nth order system is 

«) 2 (to = 0) = 1+ 2 ^ (33) 
nfc 



and the to — > 00 limit is 



^>^) = CT (34) 

The n dependence of the stability criterion is again indicative of the cut-off power-law distribution of the moments. 

It should be noted that while the natural oscillation frequency of the spring, yjk/m varies between two values, the 
above model is not at all equivalent to Kubo's model [37 of an oscillator with random frequency. In fact, Kubo in 
his paper introduced two models, one of which is described by the Kubo equation 

<j) = co{t) (35) 




FIG. 6. Time trace for a single run of the harmonic oscillator with random damping, showing the intermittent behavior. Here, 
k = 1, A = 0.2, m = 1, 7 = 0.025, a = 3. 
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where ui(t) — ujq + r), and rj is a fluctuating noise term. Here, there is only a single degree of freedom, the angular 
one, and so clearly there can be no instability in the amplitude of oscillation. Rather, there is simply a decay in time 
of the phase correlation. Kubo also introduced a model with two degrees of freedom 

x = -u){t)y 

y = cu(t)x (36) 

where w is a noisy frequency centered on ujq. Again here, there is no amplitude instability; rather the amplitude is 

exactly neutrally stable, as can be seen by the fact that x + iy can be exactly solved for as x + iy = Ae / u ^ ^ dt 
for some complex constant A. The lack of an instability can also be demonstrated for the dichotomous noise case by 
computing the eigenvalues of the Fokker-Planck operator for the moment system (x 2 )±, (xy)±, (y 2 )±, and verifying 
that they are all negative, except for an exact zero model corresponding to changes in the amplitude. Thus, in this 
case, without damping, the system is exactly marginal, whereas the random stiffness model is always unstable in the 
absence of damping. 



C. Random Mass 



Again, the fourth-order system is unstable for a wider range of parameters, and the sixth-order system for a yet 
wider range. For example, 



(OV->») = ^ (37) 



The instability boundaries are shown in Fig. [3] 



IV. CONCLUSIONS 



We have examined the dependence of the instability threshold of the harmonic oscillator with multiplicative noise 
on the mass of the Brownian particle. Even in the white noise limit, the dependence is very different for the cases 
of random damping as opposed to random stiffness. In the former case, the threshold increases as the square-root of 
the mass whereas in the latter case it is independent. The situation is even more intricate when dichotomous noise 
in consider. Then, the instability threshold is monotonically increasing for the case of random damping. For random 
stiffness, the dependence is non-monotonic for small damping, first rising for small mass and then decreasing, while 
it is monotonically decreasing for large damping. For the case of random dichotomous variation of the mass, the 
dependence is always monotonically decreasing. The stability criteria for the various cases were presented explicitly 
and their asymptotic limits discussed. 

The distribution function for the energy was computed numerically and was seen to exhibit a power-law cut-off at a 
value that increases exponentially in time. This is connected to the fact that the system has no nonlinear saturation, 
and so time serves as the only cutoff. It is also related to the fact that the instabilities of the asymptotic power-law 
distributions of < E >, < E 2 > and < E 3 >, etc. occur at smaller and smaller noise thresholds. An analytic 
calculation of the long-time asymptotic behavior of the distribution is an interesting open problem. 
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